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Abstract 



We apply a recently proposed Green Function Monte Carlo to the study of 
Hamiltonian lattice gauge theories. This class of algorithms computes quan- 
tum vacuum expectation values by averaging over a set of suitable weighted 
random walkers. By means of a procedure called Stochastic Reconfiguration 
the long standing problem of keeping fixed the walker population without 
a priori knowledge on the ground state is completely solved. In the U(l)2 
model, which we choose as our theoretical laboratory, we evaluate the mean 
plaquette and the vacuum energy per plaquette. We find good agreement 
with previous works using model dependent guiding functions for the random 
walkers. 

PACS numbers: 11.15.Ha, 12.38.Gc 
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I. INTRODUCTION 



Lattice regularization of quantum field theories is a powerful technique to compute non 
perturbative physical quantities. Its application to quantum chromodynamics has led to 
accurate predictions of the spectrum of quark and gluon states as well as of other phenomena 
like finite temperature phase transitions. 

The historical development of this active field of research developed along two main 
streams: the Lagrangian approach proposed by K. Wilson in 1974 [|1[] and the Hamiltonian 
formulation derived by J. Kogut in 1975 @. The two approaches are equivalent in the 
continuum limit but, for the purpose of analytical or numerical investigations, they offer 
quite different advantages. 

The Lagrangian approach exploits Feynman's old idea of computing quantum partition 
functions as sums over classical histories. Vacuum expectation values are obtained by eval- 
uating suitable statistical averages over long time trajectories. The basic object in this 
approach is the classical trajectory whose full temporal evolution must be retained. Space 
and time are treated in a symmetric way and a field configuration is defined over a space 
time discrete lattice. The identification of a d dimensional Euclidean quantum theory with 
a d + 1 dimensional statistical model is full. We remark that the Lagrangian approach is the 
current technique for the study of QCD mass spectra and the most precise studies suggest 
the advantage of working on anisotropic lattices with different treatment of the spatial and 
temporal discretisation. In particular, the use of spatially coarse, temporally fine lattices 
has greatly increased the efficiency of the numerical simulations [[|. 

These considerations lead to a renewed interest in the Hamiltonian formulation where 
only space is discretised and time remains continuous. This is a very natural approach in the 
study of low energy physics. Given the Hamiltonian H of a quantum many body model, its 
ground state is projected out by the application of the evolution operator U{t) = exp(— Iff), 
with t — > +oo, to any state with the same quantum numbers of the vacuum. Any accurate 
representation of the asymptotic behaviour of U (t) provides access to the ground state and 
to the low lying excitations. These representations perform clever repeated applications of 
H, in analytical or stochastic way, to specific quantum states that replace in this approach 
the role played by classical trajectories in the Lagrangian one. 

One of the first striking examples of this strategy can be found in H where the authors 
obtain good scaling measures of the mass gap in several non trivial models by analytical 
Lanczos diagonalization of the Hamiltonian. Many other analytical techniques have been 
developed based on similar ideas like, for instance, resummation of U(t) expansions or 
variational approximations to the spectrum. An updated list can be found in ||. 

On the numerical side, powerful Monte Carlo algorithms exist for the solution of quantum 
many body problems jTj . In this paper we are mainly concerned with those belonging to the 
so called class of Green Function Monte Carlo (GFMC). 

They may be regarded as the lattice version of the Feynman-Kac representation || for 
matrix elements of U(t) that are computed by averaging over ensembles of suitable random 
walkers. The dynamics of the walkers is determined by the kinetic part of the Hamiltonian 
that is its off diagonal matrix elements in the basis of walker states. The potential, the 
diagonal matrix elements of H, enters in the definition of a path-dependent weight which 
the walkers carry to represent their relative importance. In such an approach, the problem 
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is that the weights exponentially explode or vanish as the walkers diffuse and after a short 
time the algorithm becomes completely unfeasible. 

To solve this problem there are two standard classes of solutions: (a) the introduction 
of a guiding function for the random walkers, (b) a branching mechanism for the walkers 
population. Within the former method (see p|-|lTf for applications to lattice gauge theories), 
the measure in the path space is deformed and the simulation generates guided random walks 
according to a guiding function which is inspired by exactly known properties of the ground 
state (typically, its weak and strong coupling approximations). The disadvantage of the 
method is that it requires tuning of the functional form of the guiding function. This can 
only be achieved by a somewhat accurate knowledge of the vacuum structure. In fact, perfect 
guidance is equivalent to the full problem of determining the vacuum wave function. 

In the latter solution, branching is introduced among the walkers to damp their weight 
variance. Walkers with low weights are deleted and relevant walkers are replicated. The 
problem in this case is a rather involved management of the variable size walker population. 



In a recent paper [|T^], a simple procedure called Stochastic Reconfiguration, has been 
successfully applied to overcome this problem and perform simulations with a fixed number 
of walkers without introducing any guiding function, which can however be exploited if 
available. 

The aim of this paper is to investigate the application of this new technique to lattice 
gauge theories. In particular, we apply the method to U(l) 2 lattice gauge theory as a simple 
theoretical laboratory where it is easy to present the main ideas and discuss the systematical 
errors introduced by the algorithm as well as their control. 

The plan of the paper is the following. In Sec. [Tl] we review the GFMC with Stochastic 
Reconfiguration. In Sec. |IT| we apply the algorithm to the U(l) 2 model. We discuss its ac- 
tual implementation, the optimization of its free parameters and the numerical simulations. 
Finally, in Sec. |IV| we summarize our results and discuss the perspectives in the study of 
realistic gauge models. 



II. REVIEW OF GFMC WITH STOCHASTIC RECONFIGURATION 

The description of GFMC with Stochastic Reconfiguration is equivalent to present the 
Feynman-Kac formula on a lattice. For this reason we begin by discussing the simple ex- 
ample of quantum mechanics in flat space. Let us consider the one dimensional Schrodinger 
hamiltonian for a unit mass point particle 

H = \ V 2 + V(q) = H + V(q) (2.1) 

and the problem of computing its ground state wave function. We define a discrete Markov 
chain as follows. Let e be the time step (not necessarily infinitesimal) associated to each 
Markov jump and let the state of the chain be specified by the position eigenvalue q. Let 
the transition function p(q' — > q") be 

p(q'^q")=K (q",q',e) (2.2) 

where Kq is the propagator of the free hamiltonian Hq 
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K (q", q>, t) = (q"\ exp -tH \q') = -L exp J q " J? (2.3) 

An ensemble of walkers can be described at step n by its probability density P n (q)- It evolves 
according to 

P n+ i{q")= I dq' P n (q')K (q",q',e) (2.4) 
jR 

We can identify P n (q) with the wave function of the abstract state \P n ) satisfying 

Pn(q) = (q\P n ) (2.5) 
and therefore, as is well known, we obtain 

|P n ) = e-"^°|P > (2.6) 

for any finite e. 

To make a similar construction with H replaced by H we need an extension of the 
formalism. We consider a Markov chain where the state is extended from q to the pair (g, u) 
where a; is a real weight whose dynamics we shall describe in a moment. The transition 
kernel is assigned as follows 

P {q'u' -> q"u") = K (q", q', e)5 (u" - u'e'^) (2.7) 

with the correct normalization 

r roc 

/ dq" / dco" piq'u' -> q"uj") = 1 (2.8) 

JR JO 

In other words, at each discrete step, q diffuses making a random step with variance ({Sq) 2 ) = 
e as before and the weight u is multiplied by the exponential factor exp(— eV(q)). 

The probability distribution P n (q, uj) at the n-th iteration evolves according to the equa- 
tion 

P n+1 (q", u") = dq' dJ P n (q', J) p(q' , uj' - q" , uj") (2.9) 

JR JO 

To recover a wave function evolving according to exp(— tH) we must average over the weights 
and introduce the function 

roo 

ij n {q) = / duuP(q,uj) (2.10) 
Jo 

that satisfies 

iJ n+1 (q")= I dq^ n {q')e-^')K {q",q',e) (2.11) 

JR 

If we now write ip n in terms of the associated basis independent abstract state 

MQ) = (QWn) (2.12) 

we easily see that 
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|Vwi> = e- eHo e-^V„) (2.13) 
and therefore, by Trotter-Suzuki formula, in the limit e — > we obtain 



JmV t / e (g) = (?|e- w |Vo) (2-14) 

that is the Schrodinger evolution is completely reproduced. 

The meaning of the previous manipulations is that we can numerically compute the 
Schrodinger semigroup by averaging over walkers which diffuse according to the kernel K 
and which carry an additional weight uo. The weight takes care of the potential and is 
the actual relative weight of the walkers. As discussed in the Introduction, in the actual 
implementation of this method one has to face the rapidly increasing variance of the weights 
of an ensemble of walkers as the evolution time goes by. 

Stochastic Reconfiguration is a solution to this problem based on the observation that 
there are many P(q,u) giving rise to the same physical wave function ip(q). In particular 

POO fOC 

/ duo ujP(q,uj) = I duj ujP(q,oj) (2-15) 
J o Jo 



where 



P(q, u) = 5{uj - 1) / drf to' P{q, uj') (2.16) 
Jo 



and P(q,u) and P(q,u) give rise to the same wave function ip{q). The difference between 
them is hidden in the weight statistics that is not observable. The advantage of choosing 
the representative P(q,u) among all the P(q,u>) associated to a given ip(q) is that it has 
exactly zero weight variance by construction since all the weights are fixed to unity. 

The actual implementation of Eq. ( |2.15|) in a numerical algorithm is straightforward. 
Let us consider an ensemble of K walkers (characterized by their position eigenvalue and 
weight) 

£ = {(qk,u k )}k=i,...,K (2.17) 

In the K —>■ oo limit we can associate to the ensemble S a unique well defined distribution 
function Pg(q,u). We now build a new ensemble £ with K walkers and with the property 
that when K — > oo we have 

Pz(q,u) = P £ (q,u) (2.18) 

The new ensemble is simply built by extracting walkers from £ and assigning them a unit 
weight. The walkers with position q^ are extracted with probability proportional to Uk- In 
other words we extract the new K walkers from the set of old states (the values {qu}) with 
probabilities 

Pk = ^~ (2-19) 

Repetitions may occur during this multiple extraction; they change the relative frequency 
of the old low and high weight walkers. To see this, let us consider an ensemble of walkers 
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{(qi, u)i)}i=i t ... t K where for simplicity we assume that all are distinct. We want to compute 
the statistics of the fractions = rii/K of walkers with state in the reconfigured ensemble. 
Since each new walker is extracted independently, the probability of building a new ensemble 
with rij walkers in the state is given by the multinomial distribution 

if' 

p(m, . . . , n K ) = — Pi 1 ■ • • Pk (2-20) 

nil ■ ■ -n K \ 

The average number of walkers with state qi is 

(nt)= E , ' ] P T---Pk k n t = K Pl (2.21) 

niH \-tlk=K 1 

and the mean product UiUj with i ^ j 

(mrij) = K(K - l)p iPj (2.22) 

We conclude that 

(ui) = pi, (viVj) - {vi){vj) = - j^PiPj (2-23) 

Hence, in the reconfigured ensemble the state g« appears with a frequency which is precisely 
Pi ~ uji. Finite size correlations between the new walkers are present vanishing with K 
as K~ x . They induce systematic errors in the measurements that must be eliminated by 
extrapolation to the limit K — ► oo. 

We conclude this Section with a comment about Importance Sampling. The aim of 
this paper is to show that it is possible to perform simulations with a fixed size population 
of walkers without any a priori knowledge of the vacuum state |0). However, it must be 
emphasized that whenever a trial wave function for |0) is available, it can be easily included in 
the algorithm by a unitary tranformation of H as discussed in [I3|. This is straightforward 



in the Hamiltonian formulation and therefore analytical approximations (e.g. variational 
calculations) can be exploited to accelerate convergence. 



A. Generalization to other models 

From the above discussion, it should be clear that the extension of GFMC with Stochas- 
tic Reconfiguration to a more general Hamiltonian is possible only when some structural 
conditions are true. To be definite we require the existance of a complete explicit basis {\s)} 
such that H can be written 

H = T + V (2.24) 

where: (a) V is real and has vanishing off diagonal matrix elements, (b) T is the generator 
of a Markov process. Condition (b) is simply the statement that the evolution operator 
U(t) = exp(— t T) exists for t > and the kernel 

K{s\s',t) = (s"|exp(-t T)\s') (2.25) 
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is positive and normalized 

5XaV,f) = l (2.26) 

s" 

In more physical terms, we require to be able to write the Hamiltonian as a potential term 
plus a good kinetic term like, on a generic manifold, the Laplace-Beltrami operator which is 
associated to random walks on the manifold. Interesting models which belong to this class 
in the Hamiltonian formulation are the non linear O(N) a model and the SU(N) pure gauge 
theories. 



B. Observables Measurements 

In this Section we discuss how observables can be measured. The ground state energy 
is the simplest observable to be computed. To measure it, we take two arbitrary states \x) 
and with non zero overlap with the ground state |0) and write 

E ° ~ & ( X \e-*»\1>) (2 ' 27) 

The state \ip) describes the statistics of the initial state of the chain. A convenient choice, 
although not always optimal, is to take for \ip) one of the walker states and to start accord- 
ingly the Markov chain always from the same state. About \x), a good choice is to take the 
zero momentum state which is annihilated by the kinetic (off-diagonal) part of H . In other 
word, this is the generalization (on non flat manifold) of the constant wave function. In this 
case we drop the kinetic part of H and obtain 

E = lim 7, 'V ( 2 - 28 ) 

t^+oo (x\ e WI 

where V is the diagonal part of H. From the point of view of the algorithm, the projection 
over \x) is nothing but the recipe of summing over all the walkers with no additional final 
state weight. 

These prescriptions may be refined in model dependent ways, but we shall see that, at 
least in the model under consideration, they work properly. 

Concerning other operators, let us analyze in some details the measure of vacuum expec- 
tation values of operators Q which are diagonal in the chosen basis (the position eigenstates 
in the Schrodinger example). They can be computed as the limit 

= lim 0(r) (2.29) 

(0|0) T-++OC K ' K ' 

where 

«M = Jfe liffSf = <0|Q|0> + °< e " < *~* ) > (230) 

where E\ is the energy of the first excited state |1) with non vanishing matrix element 
(1|Q|0). In the above expression, the limit over t is performed automatically by the running 
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of the Markov chain. It can be considered reached as soon as equilibrium in the chain is 
achieved. The second limit requires some care and to evaluate it we allow the walkers to 
diffuse for an additional time r after the measurement of the observable. In the following 
Sections we shall check the rapid exponential convergence of Q(t) with increasing r. Of 
course, from the r dependence of Q(t) the (finite volume) mass gap Ei — E can be extracted. 



III. APPLICATION TO THE U(l) 2 MODEL 

In this Section we make the previous discussion more detailed by examining a specific 
example, the U(Y)2 lattice gauge model. Its Hamiltonian is given by 



L 

P =i 



1 d 2 



2/3 m 



(3(1 — cos i 



(3.1) 



where the link phases 9i. are defined in Fig. 1 and the gauge invariant plaquette phase <j) p is 

0p = &l,p + 02,p+l — &3,p — &2,p (3.2) 

We assume periodic boundary conditions. This lattice model has no continuum limit be- 
cause its correlation length in lattice units remains finite for all values of the coupling f3. 
Nonetheless, it shares many features with the more realistic models in higher dimensions 
and serves to illustrate the method in an easy case. We are interested in the numerical 
calculation of the ground state energy per plaquette E Q /N and of the mean plaquette 

(U p ) = (O|cos0 p |O) (3.3) 

For these two observables, the first terms in the weak and strong coupling expansions are 



known JT0| 



T 



t 0.9833 - 0.1209/T 1 + C(/?" 2 ) 



(N > 5) 



(3.4) 



(U) 



¥ 2 ~ ioP 6 + 0((3 l 



0.4917/?" 1 + 0{(3~ 



(3.5) 



A. Optimization of the Algorithm 

Let us discuss in this Section the tuning of the algorithm. Its free parameters are 

e, K, r, r (3.6) 

where e is the time step in the application of exp(— eH), K is the size of the walkers ensemble, 
r is the number of Markov chain steps between two reconfigurations and r is the time which 
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we discussed in Sec. |I B| . The best choice of r of course depends on the chosen observable. 
In principle, the same holds also for the first three parameters, but for simplicity, we shall 
discuss their optimization on general grounds independently on Q. 

The parameter e sets the scale of the elementary fluctuations of the link phases 9. Since 
((56) 2 ) = e we require y/e ~ 27r/ng with large rig to obtain a good approximation of the 
continuum diffusion. In our simulations we choose the conservative value rig = 50 and 
therefore e = 0.015. 

The value of K cannot be fixed. Instead, it must be varied and an extrapolation to 
K = oo is needed. The following functional form for all K dependent quantities 

Q(K)=Q(^) + ^ + ... (3.7) 

turns out to be enough general to reproduce quite well the K dependence. In our simulation 
we extract the three constants Q(oo), c and a from Q(K) at K = 10, 100, 1000, 5000. 

The parameter r controls the frequency of the reconfiguration process. A small r is 
useless and expensive. Actually, reconfiguration must be performed only when the walkers 
begin to show a significantly large weight variance. In that case, reconfiguration is effective. 
If the variance is small, then reconfiguration reshuffles the ensemble without useful effects. 
On the other hand, if r is taken too large, reconfiguration with a finite K will destroy the 
information contained in the ensemble and the systematic error at fixed K will be large. Our 
proposal for a choice of an optimal r is to fix it by looking at the integrated autocorrelation 
time of a relevant observable, say the vacuum energy. We set r = 1 and run the algorithm: 
the energy measurements {Ei} show an exponential decorrelation 

IS - E 2 ~ Ae~ ln ~ ml/T (3.8) 

where as usual bars denote the average over the measurements. We then set r = r. This 
procedure has two advantages: (a) the energy measurements taken after each reconfiguration 
are decorrelated, (b) the autocorrelation r is independent on K. Statement (a) holds by 
construction. Property (b) follows from the fact that r is an intrinsic feature of the chaotic 
evolution of each single walker and has no reason to be K dependent as we have checked 
explicitly. This is a useful property for the K — > oo extrapolation. 

In our simulations we consider the U(l) 2 model at (3 — 0.5, 1.0, 1.5, 2.0, 2.5 with L = 8 
to reproduce the results of | JL0[| . In this case we find that the optimal r(/3) is roughly 

r(0.5) = 40, r(1.0) = 30, r(1.5) = 25, r(2.0) = 25, r(2.5) = 20 (3.9) 

We now present our results for the ground state energy per plaquette and mean plaquette. 
For each data point we performed 10 5 Markov chain steps. 



B. Measure of Eq 

We have computed the vacuum energy per plaquette at the four values K = 10, 100, 1000 
and 5000 on a L = 8 spatial lattice. A fit with the functional form in Eq. ( |3.7| ) determines 
the /3 dependent exponent a(f3). In Fig. (|[) we show the linear fit of the numerical measures 
plotted as functions of 1/K a ^\ The results of the fit are collected in Tab. dD and compared 



9 



with the weak and strong coupling expansions. The results are quite good, especially because 
the strong coupling series is only known up to the next to leading term. The CPU time 
required for each point is about one hour on a Pentium 200 processor. Moreover, we stress 
that to obtain these results we only needed a very simple tuning of the algorithm and no 
additional knowledge of the ground state properties. As can be seen from the figure, the 
extrapolation is very near the value obtained with the largest K used. In the present model, 
it is fairly easy to obtain analytical approximations for the ground state energy per plaquette. 
The hamiltonian written in terms of the gauge invariant phases (f) p is 




N 

and the independent plaquette approximation gives 

j-i(Mathieu) -, 



(3.10) 



/?+7^ a o(-/n (3.H) 



L r 2/3 

where ao(q) is the lowest characteristic value for the even solutions of the Mathieu equation 

y"(x) + (a (q) - 2qcos{2x))y = 0. (3.12) 

A more practical approximation can be obtained by a variational calculation based on the 
following simple independent plaquette Ansatz 



V>q(0i, ■■■An) = exp I A ]T (■<>*'">,, j • (3.13) 

In this case we obtain 




E^ r \p) _ 1 . 



(A-/3 



2 /i(2A) 
j / (2A) 



(3.14) 



(where I n is the n-th modified Bessel function). In Tab. (|I|), columns Mathieu and var show 
the numerical values of these approximations. 

As one can see, they appear to reproduce well the Monte Carlo data over the whole 
region of the coupling. This holds true especially for the estimate computed in terms of 
Mathieu functions. However, we stress again that the Monte Carlo extrapolated data have 
no systematic errors and have been obtained with the minimum a priori knowledge. As 
stated above, any analytical information on the ground state, like a variational ground state 
wave function, can be used for the acceleration of the Monte Carlo simulation by including 
in the algorithm an importance sampling step. 



C. Measure of (U p ) 

The results obtained for the average plaquette are quite similar to the previous ones. 
Again, we show in Fig. @ the linear extrapolation after the determination of a(/3). In 
Tab. (0) we see that the Monte Carlo data is well matching the analytical series except 
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in the region around (3 = 1.0. However, the weak coupling series is not to be trusted at 
(3 = 1.0 since its plot shows a steep variation between (3 = 1.0 and (3 = 1.5 signalling the 
need for an additional term. The point at (3 = 1.5 deviated by about 3% with respect to the 
strong coupling series. This can simply explained by assuming that it is a too small value 
for strong coupling to apply and moreover the next points at larger (3 match better, below 
the percent level. In Fig. we show the dependence of (U(t)) on the time r to check 
the exponential convergence to the vacuum expectation value. Again, we stress that these 
numbers are obtained without a priori information and over the whole coupling variation. 

To summarize our results we collect in Fig. (^) a graphical comparison of Monte Carlo 
data with the asymptotic (3 — > and (3 — > oo expansions. 



D. Computational Cost 

In this Section we address the problem of estimating the computational cost of the 
algorithm. A realistic computation must follow two steps: 

a) extrapolation K — > oo at fixed coupling f3: 

b) continuum limit — > oo where £ is the correlation length in lattice units (we do 
not discuss finite volume effects and their control). 

In performing step (a) we need the computational cost as a function of K that is the CPU 
time required to achieve a give statistical error e sta t in the evaluation of an observable Q. 
With N Markov chain steps we have 

tstat — — 1= - (3.15) 

IN/r((3) 



where ctq(/3, K) is the standard deviation of Q measurements and N/r(f3) is the number of 
independent measurements. The dependence of oq on K is simple. Due to a self average 
effect we expect ctq(P,K) = aq(P)K~ 1 ^ 2 . Our explicit numerical simulations confirm this 
scaling law. The CPU time is roughly proportional to NK since the reconfiguration process 
turns out to be only a small fraction of the computational cost. Therefore 

e sta t ~ , = (3.16) 

^/T CPU /r((3) 

and K cancels. The increase in CPU time associated to the management of the ensemble is 
compensated by the larger statistics. In other words, the cost of stochastic reconfiguration 
is just the cost of the fit described by Eq. ( |3.7| ): it is proportional to p the number of k 
values used to extract Q(oo) from Q(K). 

The cost of step (b) is expected to be much more model dependent. As the continuum 
limit is approached, the quantities which vary with (3 are &(/3), r(/3) and a((3). There is no 
mechanism in this algorithm to prevent slowing due to increasing r{(3) when the correlation 
length increases and the behaviour of the algorithm must be studied in the case of realistic 
models with diverging correlation length. What we observe in the considered model is a 
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very mild j3 dependence of the exponent a(/3) suggesting that Eq. ( |3.7|) may hold without 
the need for a huge K. Moreover, an advantage of the Hamiltonian formulation is that 
this slowing could be eventually reduced by importance sampling as discussed in the end of 
Sec. H; in fact, preliminary investigations show that the use of a good trial wave function for 
the ground state can strongly reduce and allow to perform simulations with reasonable 

errors using relatively small values of K. 



IV. CONCLUSIONS AND PERSPECTIVES 

In this paper we applied a recently proposed many body Monte Carlo algorithm of 
Green Function type to the study of lattice gauge models in Hamiltonian form. This kind 
of algorithms compute the vacuum wave function and many related quantities by averaging 
over a set of suitable weighted random walkers. The method that we discussed solves the 
annoying long standing problem of fixing the walkers number without exploiting a priori 
informations on the ground state structure and is thus quite general. The trade off is 
the introduction of a systematic error, but we showed by explicit simulations in the U{1)2 
model that it can be kept completely under control. The algorithm can be applied to more 
interesting models like the non linear O(N) a model or SU(N) pure gauge on which work 
is in progress |14| . The perspectives of this works seem interesting expecially in view of the 



recently proposed improved lattice Hamiltonians fl5j for pure gauge models where we plan 
to investigate the algorithm efficiency. 
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FIGURES 

FIG. 1. U(l) N-chain geometry. Definition of the link angles. 
FIG. 2. Extrapolation K — > oo for the vacuum energy per plaquette. 
FIG. 3. Extrapolation K — ► oo for the mean plaquette. 
FIG. 4. Large r limit of (U(t)). 
FIG. 5. Graphical comparison of Monte Carlo data with strong and weak coupling series. 
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TABLES 



TABLE I. Eq/L compared with weak and strong coupling approximations. The last column 
shows the relative deviation between the Monte Carlo (MC) data and the weak or strong expansions. 



(3 


weak 


MC 


Mathieu 


var 


strong 


% A 


0.5 


0.4689 


0.4694(2) 


0.4690 


0.4690 




0.1 


1.0 


0.7732 


0.7736(2) 


0.7724 


0.7746 




0.05 


1.5 




0.8907(2) 


0.8909 


0.9005 


0.9027 


1 


2.0 




0.9292(2) 


0.9299 


0.9435 


f\ noon 

0.9229 


0.7 


2.5 




0.9524(2) 


0.9466 


0.9594 


0.9349 


2 




TABLE II. (U) compared with weak and strong 


; coupling 


approximations . 




P 


weak 




MC 




strong 


% A 


0.5 


0.1236 




0.1229(2) 






0.6 


1.0 


0.4073 




0.4262(2) 






5 


1.5 






0.6513(2) 




0.6722 


3 


2.0 






0.7590(2) 




0.7542 


0.6 


2.5 






0.8074(2) 




0.8033 


0.5 
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